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ABSTRACT 

Using the theoretical model of E. A. Pashitskii et al. We have 
studied the various semi empirical potential for interaction of 
helium atoms in real space. We have also theoretical evaluated the 
Bogolyubov quasi particle spectrum in Bose liquid and compared 
our results with inelastic neutron scattering experiments. Our 
theoretical results are in good match with the experimental data. 
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1. INTRODUCTION 

Although there are large numbers of 
calculation 1 6 for the spectrum of elementary 
excitation in this super fluid (SF) 4 He Bose 
liquid but there remains some problem in 
this direction. There is an excellent 
agreement will experimental data in the 
region of the roton minimum obtained by 
Monte carlo method making use of the 
shadow wave function 1 and by the 
correlation basis function 2 . These calculation 
also employed the modern inter atomic 
potentials 35 for 4 He. There are some 
calculation using microscopic perturbation 
theory 6 8 for long wave phonon part of the 
spectrum E(p) - Cip where Ci is the speed 
of the first ( hydrodynamic) sound in liquid 



4 He. These calculation faces some principal 
difficulties become the no renormalized 
perturbation theory gives rise to infrared 
divergence and are non-analytical 9 " 12 at p — > 
0 and e — > 0. These difficulties have been 
removed by the application of combined 
variable techniques 13 . These variables 
reduces to hydrodynamics variable of 
macroscopic quantum hydrodynamics 14 in 
the long-wavelength limit (p — > 0). On the 
other hand in the short wavelength domain 
they correspond to the bosonic quasi particle 
creation and annihilation operators. 

Neutron inelastic data 15 ' 1 and results 
in quantum evaporation of 4 He atoms 17 so 
that the maximal density p 0 of the single 
particle BEC (Bose -Einstein Condensate ) at 
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T < Tk does not exceed of the total density 
p of the liquid 4 He, whereas the density of 
the SF component p s — > p and T — > 0. This 
indicates that there is strong interaction 
between 4 He atoms and the quantum 
structure of the super fluid condenset in He- 
ll carry an excess density (p s - p 0 ) » p 0 . 
This requires much more investigations. 

In this paper, we have studied 
various semi empirical potentials to study 
helium-helium interaction. These potentials 
involve strong repulsion at small distances 
and weak Van-der Waals attraction at large 
distances. We have also computed the 
Bogolyubov spectrum of a dilute quasi ideal 
Bose gas and compared our evaluated results 
with inelastic neutron scattering data. Our 
evaluated results are in good match with the 
experimental data. 

2. MATHEMATICAL FORMULA USE 
IN THIS STUDY 

In order to describe the interaction 
of He atoms in real space one uses various 
semi empirical potentials. These potentials 
involve strong repulsion at small distances 
and week Van-der Waals attraction at large 
distances. However, most of those potentials 
are characterized by strong divergence at r 
— > 0. For example Lennard-Jones 
Potential 18 . 





12 




\r ) 







e and a are potential parameter r>r c . 
If one takes the Fourier transform of the 
potentials 

V(p) = jd 3 rU(r)e 2pr 



= — [rU(r)Sin(pr)dr 

P ° (2) 

One cannot use this potential for the 
description of pair interaction in nonnative 
space as it is infinite, diverging as the lower 
limit. 

The other potential which is used for 
the calculation of inter atomic interaction 
and of possible bond state is Aziz 
potential 19 . 

(/ A (r)=Ae Xt X^-A 2 )-expK(^/r)-l) 2 ]XQ^- 2 "r< ?;) 

J5=0 

= Aexp(-«r -/3r 2 ) - ^ C 2k+6 r~ 2k ~ 6 r > r 0 
Where, 

A = 1.8443101 x 10 5 K 
a = 10.43329537 A 1 
P = 2.27965105 A 2 
C 6 = 1.36745214 Kx A 6 
C 8 = 0.42123807 K x A 8 
Cio = 0.17473318 K x A" 10 

Such potentials remain finite at r = 0 
due to the non analytic exponential 
dependence on r which suppresses any 
power divergence at r — > 0. The potential (3) 
is convenient for calculations in real space 
making use of J as trow-like wave functions. 
But employing its Fourier component for 
solving nonlinear integral equations is 
technically difficult. The model potentials 
which describes interaction of helium atoms 
in real space fact to take into account that 
the distance less than the quantum radius of 
the helium electron shell r 0 = 1.22 A the 
Coulomb repulsion between the nuclei sets 
in. 
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Piashitskii et al. 2Q suggested a 
potential which diverges as r 1 at r — > 0 



42 2 f 
U(r) = (l + //r)exp 



V 



L] 
a) 



r<r 



momentum transfer p. the same oscillating 
Fourier component is characteristics of 
smooth potential V(r) in the form of 
Linhardt type function 22 having and infinite 
negative derivatives at the inflection point 
r = a. 





12 


[-] 


\r j 







r>r c (4) lK 2 



From the conditions of continuity of 
the potential U(r) and its first derivatives at 
r = r c , one determine the value of the 
parameters r c = 2.38A and |i = 229 A 1 for 
a = r 0 / 2 = 0.61 A, e = 10.8 K and a = 2.642 
A. 

In many body problem, the quantum 
effects play an essential role and in that case 
the two-body potential is irrelevant. For that 
one should have a model potential with a 
simpler analytic expression. There is a 
model potential in the form of a Fermi type 
function in real space 



U F (r)=V 0 



exp 



V-q 2 ^ 



+ 1 



(5) 



Which at b = 0 degenerates into a 
state function at 0 (a-r) with a finite height 
U 0 at r < a. this model corresponds to a 
model of 'soft' sphere. Its Fourier 
component is expressed in terms of first 
order spherical Bessel function 21 

V(p)=V 0 ^, fl (x) = ®»W-f°«*r l 
pa x 

(6) 

Where V 0 = 3V(0) = 471 U„ a 3 . It is an 
oscillating sign changing function of 



1 + - 



2r 



a + r 



(7) 

Formally this problem is an inverse 
to one of periodic oscillations of spin density 
in real space of the exchange interacting 
spins of electrons in metal. This is popularly 
known as RKKY (Ruderman-Kittel- 
Kaisuya-Yosida) oscillation 23 . The same 
behavior is characteristics of the Fourier 
components of more realistic potentials 
diverging faster than r 1 at r = 0. It posses 
inflection points in the radial dependence at 
r = r c . The amplitude of the oscillation of the 
Fourier component of the Fermi type 
potential (5) at b ^ 0 is falling off 
exponentially with the increase of the 
parameter b due to the decreasing absolute 
value of the negative derivative at the 
inflection point. 

If our situation the oscillation 
potentials (6) into the Bogolyubov of a 
dilute quasi-ideal Bose-gas 



24.25 



E B (P) = 



2m 



+ 2nV(p) 



(8) 



3. DISCUSSION OF RESULTS 



In this paper in order to describe the 
interaction of He atoms in real space, 



Journal of Pure Applied and Industrial Physics Vol.2, Issue 3, 1 July, 2012, Pages (142-285) 



207 Md. Shabbir, et al, J. Pure Appl. & Ind. Phys. Vol.2 (3), 204-210 (2012) 



various empirical potentials are used which 
involve strong repulsion at small distances 
and week Van der waals attraction at large 
distances. Most of the potentials are 
characterized by a strong divergence at r— >0. 
Among the various potential discussed, there 
is a repulsive potential in the frame work of 
'soft spheres model' the Fourier component 
v(p) is an oscillating sign - changing 
function of momentum transfer p. It appears 
that in a certain region of momentum space 
at p ^ 0 there is an effective attraction 
between bosons v(p) < 0. This has nothing to 
do with Van der waals forces and has a 
quantum mechanical diffraction nature. This 
attraction gets substantially enhanced due to 
the multi particle effects of renormalization 
(screening) of the initial interaction 28 30 . 

We have theoretically compared the 
Bogolyubov quasi particle spectrum 
equation (8), with the oscillating Fourier 
component [v(p) = V 0 Sin (pa) / pa] and also 
with experimental spectrum 26 ' 27 . The results 
are in close appearance with the hard sphere 
potential (equation 6). The results are shown 
in table T 5 . In table Ti we have shown the 
evaluated radial dependent of potential (4) 
and Aziz potential (5), in table T 2 the two 
forms of finite potential (5) has been shown 
in real space (r/a) keeping b = 0 [I - form] 
and b = 0.5a [II - form]. In Table T 3 . We 
have given the (Lindhard function potential 
[U L (r) / Uo] shown in equation (7) in real 
space (r/a). In Table T 4 , we have shown the 
momentum dependent of the Fourier 
component of the Fermi type potential (5) 
for b = 0 and potential (4). 

From these studies it looks that the 
soft sphere model repulsive potential which 
gives very good agreement between the 
theoretical quasi particle spectrum E B (p) and 



the experimental spectrum of elementary 
excitation in quantum Bose liquid 4 He, is 
much smaller than the value of the Aziz-type 
potential at r— >0. It is a result of strong 
quantum diffraction effects in Bose liquids, 
because the average distance between 
particles is equal or less than the de Broglie 
wavelength for bosons. 

Table T l 

An evaluated results of radial dependence of 
potential (4) and Aziz potential (3) 



r(A) 


Potential (4) 
U (r) (K) 


U A (r)(K) 


2.5 


105.2 


110.8 


2.6 


95.6 


106.5 


2.7 


84.3 


98.2 


2.8 


50.6 


68.6 


2.9 


16.7 


23.2 


3.0 


-12.2 


-16.8 


3.1 


-14.6 


-15.6 


3.2 


-10.4 


-11.8 


3.3 


-8.6 


-9.9 


3.4 


-6.5 


-7.3 


3.5 


-5.6 


-6.8 


4.0 


-4.32 


-5.8 


4.5 


-4.00 


-3.8 


5.0 


-0.08 


-1.02 
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Table T 2 

An evaluated results of the finite potential 
model (5) in real space (r/a). [I - form for b=0] 
and [II - form for b=0.5a.] 



Table T 3 

An evaluated results of the Lindhart function 
potential U L (r) / U 0 equation (7) in real space 
r/a. 



r/a 


U F (r) / U 0 for 
b = 0 


U F (r) / U 0 for 
b=0.5a 


0.0 


1.00 


1.00 


l al 


1.00 


0.998 


0.2 


1.00 


0.996 


0.3 


1.00 


0.987 


0.4 


1.00 


0.965 


L 0.5 


1.00 


0.928 


0.6 


1.00 


0.889 


0.7 


1.00 


0.793 


0.8 


1.00 


0.786 


l 09 


1.00 


0.684 


1.0 


0.00 


0.546 


1.1 


0.00 


0.253 


1.5 


0.00 


0.446 


1.6 


0.00 


0.054 





T T ff\ 1 T T 


u.u 


1 .uu 


0.1 


0.984 


0.2 


0.975 


l °- 3 


0.952 


0.4 


0.886 


0.5 


0.825 


i 0.6 


0.795 


0.7 


0.705 


0.8 


0.659 


l °- 9 


0.613 


1.0 


0.526 


1.1 


0.508 


i 1.2 


0.432 


1.3 


0.326 


1.4 


0.286 


L L5 


0.239 


1.6 


0.187 


1.8 


0.144 



Table T 4 

An evaluated results of the momentum dependence of the Fourier components of the 
Fermi type potential (5) for b = 0 and potential (4). 



P (A" 1 ) 


V(p)/V(o) 


Potential (5) 


Potential (4) 


0.0 


1.0 


1.0 


0.2 


0.869 


0.986 


0.4 


0.792 


0.924 


0.5 


0.607 


0.885 


0.6 


0.496 


0.816 


0.8 


0.288 


0.774 


1.0 


0.087 


0.625 


1.5 


-0.156 


0.463 


2.0 


-0.108 


0.136 


2.5 


-0.088 


-0.115 


3.0 


0.145 


-0.056 


3.5 


0.208 


0.069 


4.0 


0.125 


0.145 


4.5 


0.086 


0.186 


5.0 


-0.095 


0.068 
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Table T 5 

An evaluated results of Bogolyobov quasi particle spectrum (equation 8) with the 
oscillating Fourier component of the hard sphere potential equation (6). 



p (A- 1 ) 


E B (p), K 


Hard sphere 
Potential eq (6) 


Bogolyobov quasi particle 
Spectrum eq" (8) 


0.0 




U.U40 




3.628 


4.554 


0.50 


6.836 


7.397 


0.75 


10.957 


11.412 


1.00 


14.254 


15.350 


1.25 


12.868 


13.145 


1.50 


10.296 


11.254 


1.75 


8.144 


9.266 


2.00 


9.686 


10.052 


2.25 


11.234 


12.645 


2.50 


16.286 


17.128 
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